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The concept of transformation optics is extended to nonlinear electrodynamics. It is shown that 
transformation optics favors implicit constitutive relations in terms of energy densities D ■ E and 
B ■ H rather than E^ and H^. The Kerr nonlinearity is studied in detail and the transformation 
optics based engineering of self-interaction effects is discussed. As a specific example we introduce 
transformation optics applied on a self-focusing field configuration. 



o 

O 
D 

Q 

in 



O 

Oh 

q 
o 

Oh 



> 



1^ 

o 



^ 

cd 



I. INTRODUCTION 

Transformation optics (TO) [IH^ has become one of 
the most powerful tools in the design of artificial media 
(metamaterials) . TO allows the derivation of a linear 
constitutive relation directly from a desired trajectory of 
light - in a completely algebraic way. The key ingredient 
of TO is to define the trajectory in the specified medium 
as a result of a transformation applied on the trajecto- 
ries in free space. Unlike normal coordinate transfor- 
mations, under which physics is invariant, the field lines 
are considered to be attached to the coordinates during 
the transformation, thereby changing the physics and de- 
riving the constitutive law of the desired medium. De- 
spite some important differences, which are explained in 
Sect. [ll| the characteristics of transformation media re- 
main very close to those of the transformed empty space. 
Thus, transformation media are often interpreted as me- 
dia that mimic (a transformation of) empty space. This 
possibility greatly stimulated research in metamaterials: 
besides the original proposal of a cloak [H [2] (a point is 
blown up to a sphere) , the large number of proposals in- 
clude beam shifters and benders |4j, field concentrators 
[S] (a large cylinder is compressed into a smaller one), 
flat hyperlenses by means of stretching coordinates [6l 
and, more recently, a carpet cloak ■?;. Also, attempts 
have been made to generalize the concept of TO: instead 
of a Euclidean empty space a curved space can be mim- 
icked [8 , which can be used to reproduce effects known 
in general relativity |^ , and more complicated geometric 
[TU] or even direct field transformations [TT] can replace 
the simple coordinate transformation. It should be men- 
tioned that TO also stimulated research outside the field 
of electromagnetics, notably similar principles can be ap- 
plied to matter waves in quantum mechanics |12) or to 
sound waves in acoustics [131 El] ■ An extensive review of 
the mathematical principles of TO can be found in Ref. 
[15j , a recent summary of various applications has been 
presented in Ref. [TB]. 

Since transformation media are linked to empty space 
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via a coordinate transformation it follows immediately 
that all transformation media must be linear. Trans- 
formation of time in standard TO leads, in general, to 
transformation media with time dependent media param- 
eters. In these media effects such as frequency conversion 
are possible |17j . although the material is linear. Similar 
conclusions apply to the aforementioned extensions ex- 
cept for the field transforming metamaterials [TT], which 
however lack a simple geometric interpretation. More- 
over, it was shown recently that the constitutive relation 
obtained from TO is unique up to some rescaling factors 
[T8] . In this paper the idea of TO is applied to nonlinear 
electrodynamics. From the conclusions drawn above the 
basic strategy of this task is immediate: if a nonlinear 
medium is to be obtained after the transformation, one 
has to start with nonlinearities. Therefore, a nonlinear 
medium cannot be seen as a medium "mimicking free 
space" as is the case in standard TO. Still, the transfor- 
mational approach can provide a tool with which solu- 
tions in complex, nonlinear media can be derived from 
much simpler - though nonlinear - situations. The very 
generic setup of TO guarantees that the solutions of the 
Maxwell equations in the complex material follow alge- 
braically from the solutions in the simpler situation and 
- as in linear TO - this algebraic transformation has an 
immediate geometric interpretation. In this way, whole 
classes of nonlinear media can be related via TO to one 
particular, simple case (as, for example, the Kerr nonlin- 
earity). Due to symmetry properties it will be derived 
that the nonlinear interactions of these classes are de- 
scribed in terms of energy densities D ■ E and B ■ H 
rather than E'^ and H^. 

Most applications of linear TO are supposed to be in- 
dependent of the electromagnetic field configuration - as 
an example, the TO cloak 1,^2, theoretically produces an 
invisible area for any incoming electromagnetic field. All 
these applications straightforwardly extend to nonlinear 
TO. However, even simple nonlinear media such as a Kerr 
material exhibit interesting self-interaction effects, which 
depend not just on the medium parameters but also on 
the details of the electromagnetic fields propagating in 
the nonlinear medium. As is shown, such effects can also 
be engineered effectively using nonlinear TO. Of course. 



the workability of the resulting devices will be restricted 
to exactly those field configurations that exhibit the in- 
teraction effects already in the original untransformed 
medium, but the main advantages of TO (the geometri- 
cally intuitive construction and the algebraic derivation 
of the solutions) are completely retained. 



II. TRANSFORMATION OPTICS FOR LINEAR 
MEDIA 

In its standard formulation, TO is restricted to a spe- 
cial class of linear media that may be given the interpre- 
tation of mimicking the free space solutions. To facilitate 
the discussion of specific issues appearing with nonlinear 
TO, this simpler concept is introduced in this section - in 
doing so, some technical aspects important for this paper 
are highlighted. 

Since TO deals with generic, curvilinear coordinates, a 
suitable notation taking care of the covariance properties 
in such coordinates should be used. Here, we use compo- 
nent notation in conjunction with the Einstein summa- 
tion convention. The Maxwell equations in this notation 
read 



V,B* - , 



dt 



''"d.Ek - 



dt 



-D' 



J 



(1) 
(2) 



Here, di = d/dx^ is the directional and V^ is the covariant 
derivative with respect to the spatial metric metric 7ij 
and its determinant 7: 

v.A* = (a, + r^,)A' = -^a.(^A') . (3) 

Finally, e*^*^ is the Levi-Civita tensor, 



ijk 



%/7 



[ijk] ; 



)1 (ijk) are an even permutation of (123), 
— 1 (ijk) are an odd permutation of (123), 
otherwise. 

(4) 

Standard TO starts with the Maxwell equations comple- 
mented by the free space constitutive relations 



D'=eo-f'^Ej 



B'^^iof'H, 



(5) 



In free space the trajectories of light are straight lines, 
and to obtain a transformation medium a coordinate 
transformation is applied to this situation. Of course, a 
normal coordinate transformation does not affect physics, 
but in TO the electromagnetic field lines are considered 
as "attached" to the coordinates, such that the coordi- 
nate transformation does not represent an invariance but 
in fact deforms and stretches the trajectories of light. 



The resulting transformation medium is thus not physi- 
cally equivalent to the original free space situation, but 
nonetheless the complete solution of the Maxwell equa- 
tions in the medium is obtained in an algebraic way from 
the free space solution. 

To understand the technical aspects of TO, the covari- 
ance properties of all quantities involved must be under- 
stood. In index notation, these covariance properties are 
encoded in the position of the indices. D = [D^] and 
B = [B'] with upper indices are vectors (column vec- 
tors), while E = [Ei] and H = [Hi] are co- vectors (row 
vectors). By means of the scalar product a column and a 
row vector are mapped onto a scalar, e.g. D ■ E — D^Ei, 
but no similar procedure exists for two vectors or two co- 
vectors; in terms of the index notation this manifests in 
the rule that summation is only allowed over opposite in- 
dex positions. Vectors are transformed into co-vectors by 
means of the metric, Di = D^^ji, co- vectors into vectors 
by means of the inverse metric, E^ — ^^ Ej. This manip- 
ulation allows to define a scalar quantity E"^ — Ei^^^ Ej 
with an analogous expression for vectors. Since the met- 
ric in Cartesian coordinates is simply the unit matrix, no 
essential difference between vectors and co-vectors exists 
in that case, but in generic coordinates a careful distinc- 
tion of the two kinds of objects is indispensable. This 
is seen when considering coordinate transformations, the 
first step in the TO program. 

Under a transformation x — > x{x) with the Jacobian 



rpl 



dx3 



(6) 



each index transforms linearly, again following the rule of 
summation over opposite indices and the obvious require- 
ment that a summation can only be performed over in- 
dices that belong to the same coordinate system. There- 
fore [J 



&[x[x)) 
E, {x{x)) 

lij (s(a;)) 



D^{x) 



9x* 
dx^ 
dx'' dx^ 



(7) 
(8) 
(9) 



To express the effect of such a transformation on Eqs. ([I]) , 
(I2]) and ([5]) we can merely put a bar over each quantity, 
since physical processes are independent of the choice of 
coordinates. The content of the equations remains ex- 
actly the same. As mentioned above, TO makes use of 
special transformations, where the field lines of electro- 
magnetics are considered to be stitched to the coordi- 
nates. More technically this means: While the barred 



^ In this work we consider exclusively proper coordinate transfor- 
mations, det[T*j] > 0. The extension of TO to include improper 
coordinate transformations includes some technical difficulties, 
but does not yield any new types of material |19l . 



metric should be used in the constitutive law (51), the 
Maxwell equations should be written in terms of the orig- 
inal metric without a bar. This ensures that the solutions 
of the Maxwell equations indeed are written in the coor- 
dinates X with metric 7^ , but light follows straight lines 
with respect to the metric jij and thus behaves "as if" 
the metric were 7^^ . This manipulation is possible thanks 
to the fact that the metric appears in the Maxwell equa- 
tions only via the covariant derivative ([3]) and therein 
the metric 7ij can easily be changed back to the original 
metric 7^^ by rescalings [21 [12] 



D' 



V7 



zD' 



B' 






B' 



Ei = Ei 



Hi = Hi 



(10) 



space solutiorQ Besides the relative orientation of vec- 
tors and CO- vectors discussed above this applies in partic- 
ular to the stress-energy-momentum tensor |10j and the 
Poynting vector as a part of the latter. 

It should be mentioned that more rigorous, coordi- 
nate independent formulations of TO are possible. Also, 
since electrodynamics is inherently relativistic, transfor- 
mations mixing space and time can be addressed, which 
yields an extension of TO to bi-anisotropic media [3]. 
However, the interpretation of such transformations is 
less immediate in the presence of nonlinear media, thus 
we will restrict ourselves to purely spatial transforma- 
tions in this work. 



which modify the transformed vacuum relation to the 
constitutive law of TO 



III. TRANSFORMATION OPTICS FOR 
NONLINEAR MEDIA: A GENERIC ANSATZ 



D' = 



Vi 



Sof'E, 



B' ^ 



V7 



fXof'H, 



(11) 



It should be stressed again that these equations have to 

be understood as equations in the original coordinates x. 

The rescaling procedure ( 10 ) is not a symmetry of 



electrodynamics since it maps a free space solution of 
the Maxwell equations onto a solution in a certain lin- 
ear medium. Therefore, there must exist some physical 
differences between the free space and the TO solution. 
Again, the component notation proves extremely power- 
ful for clarifying some of these differences. Obviously, in 
any free space solution D is parallel to E, if the latter is 
transformed into a vector: 

{D\D',D') II iE\E^,E') = (7i'S„72»i?„73'ii;,) , 

(12) 

{D\D\D^) II {E\E\E^) = {j^'E,,f'E,,f'E,) . 

(13) 

Evidently, being parallel is a coordinate independent 
statement and thus holds in the transformed coordinate 
system x as well. The situation is different in TO: Here 
the constitutive law is written in terms of the matrix 7^^ , 
the metric in a fictuous "electromagnetic space." Still, 
co-vectors are transformed into vectors by means of the 
true metric of the physical space, which is "fij . Therefore, 

iD\D',D^) = ^e,{f=E„f^E,,f'E,) 

^(7i^^„72*^„73^^0- (14) 



As mentioned above, the only strategy for obtaining a 
nonlinear transformation medium is to choose as a start- 
ing point nonlinear medium constitutive relations instead 
of the free space relations (Is]) . In this section we analyze 
how the general ansatz of a nonlinear medium 

D^ = (e^= + x'tEk + ^t'EkEi + ..)e, , (15) 
B' = U^ + epi^fc + (f'Hr^Hi + ..)Hj , (16) 



behaves in TO. Here, e*-' = EoE/j and /i*^ = /ioMj^ are 
the - in principal arbitrary - linear permittivity and per- 
meability tensors, while the tensors Xi i xi ^l*^- and 



^1 I ^2 ^tc. parametrize the nonlinear contributions. A 
transformational approach to design new media starting 
from Eqs. ( 15 ) and ( 16 ) is straightforward. In the first 



step we apply the coordinate transformation; the nonlin- 
ear constitutive law is invariant if all medium parameters 
transform as tensors: 



e^o ^T\T^ js'^ xi 

-i j' k I 

X2 = 



i'j'k' rjii' rpj' rpk\ ^ijk 



j^ kxr (17) 



pi'j' ^T\T^' j^}^ S;[^'^' =T\T^' jT^'kCi'' (19) 



y-i j k I rrti' rrii' rrik' rril' /-ijkl 



(20) 



In the second step the rescaling rules ( 10 ) are applied. 



which do not affect the new nonlinear contributions at 



Thus, in a transformation medium D and E fail to be 
parallel in general, an observation that has important 
consequences in what follows. This observation can be 
generalized: Quantities (except the constitutive relation) 
that depend on the metric are not mapped correctly 
within TO and, thus, represent a physical difference be- 
tween the transformation medium and the original free 



^ The metric dependence of Gauss's law and its magnetic coun- 
terpart is somewhat accidental and easily can be removed in the 
so-called premetric formulation of electrodynamics [20' . In the 
present formulation this dependence can be removed by means of 
rescalings JTM since the covariant derivative just depends on the 
determinant of the metric, but not on its individual components. 



all. Thus, the result 



D' = 
B' = 






n 



'''''Hk+S,t'HkHi + . 



Ej , (21) 
H, , (22) 



is immediate. 

This result shows that TO extends straightforwardly to 
nonlinear electrodynamics. Nonetheless, the interpreta- 
tion of the result is not as simple as in standard TO, but 
it can be cast into the following statement: The transfor- 
mation medium (121) and (22 1 behaves as if it were the 



original medium ([15) written in coordinates x, but in fact 
it is the physically different medium (21)/(22| in coordi- 
nates X. Here, "behaves as" has the same meaning as in 
standard TO, namely the solution of the Maxwell equa- 
tion in the transformation medium with coordinates x is 
- up to rescalings ( 10 ) - equivalent to the solution of the 



original medium in coordinates x, but this analogy does 
not extend to any quantity that depends on the metric. 
In this way it is possible to design complicated nonlin- 
ear media starting from a simple, though still nonlinear, 
situation. 

It was mentioned above that four-dimensional space- 
time transformations instead of purely spatial coordinate 
transformations can be applied in standard TO [3] . Since 
transformations mixing space and time transform electric 
into magnetic fields and vice versa, such transformations 
are associated with bi-anisotropic media. In the pres- 
ence of nonlinearities, the mixing of electric and mag- 
netic fields affects not just the linear term, but also the 
nonlinear expansion. Thus, such transformations would 
result in nonlinear electric-magnetic mixing terms, which 
might be difficult to interpret. 



IV. EXPANSION IN ENERGY DENSITIES AND 
THE KERR NONLINEARITY 

A. The general expansion 

In many applications of nonlinear electromagnetics the 
nonlinearities entering the constitutive law, the polariza- 
tion density or derived quantities such as the refractive 
index are parametrized in terms of E'^ and H^. Instead 
of the generic law ( 15 ) and ( 16 ) this suggests the common 
ansatz 



D^^(e^^+XI'E' + X^{E 



2\2 



B' = (^I'^+^l'H^ + a{Hy 



E, 



H, 



(23) 
(24) 



where the nonlinear terms are now parametrized in terms 
of sequences of matrices Xi and f}-' (/ = 2, 4, 6, . . .). Of 
course, we could try to use the standard formulae of TO 
to derive the constitutive law of transformation media 



starting with a special case of Eqs. (23) and (24). But at 



this point, we should notice that the expansion parame- 



ters in Eqs. (23) and (24 1, 



£;2 



E.^'^E, 



B' 



Ha^'H, , (25) 



implicitly depend on the metric. From the transforma- 
tion rule 



£;2 = Ea^'Ej 



E^ = E.f^E, 



^ {E ■ T-i),(r • 7 • T'^yi{T-^T . ^)^ 



= E^ 



(26) 



we thus would find that after the transformation the re- 
sulting nonlinear medium. 



D' 



-1-j 



$,i^ 
,..^i. 



'2n2 



xl'E'+xiiE') 



-e^H^+ci{Hy 



E, 



H 



] 1 



(27) 
(28) 



is not expanded in terms of the fields in the transfor- 
mation medium E"^ — Eij''^Ej, H^ = Hij^^Hj, but in 
terms of the fields in the original medium as in Eqs. ( 23 ) 
and pll). 



This difficulty can be overcome by replacing the 
conventional expressions (|23|) and (24) with equivalent 



metric-independent expressions. To this end, we notice 
that in an arbitrary (and potentially anisotropic) mate- 
rial the energy densities are not proportional to E^ and 
H^. Instead, an invariant energy density of the electro- 
magnetic field is given by 



U 



D E + B H 



(29) 



which suggests using the electric and magnetic energy 
densities. 



Ue = D- E = D'Ei , Um = B H = B'H, 



(30) 



as expansion parameters. As shown by these equations, 
this approach automatically removes all concerns about 
metric dependence since the invariant energy densities 
are indeed independent of any spatial metric. This means 
that due to both physical and symmetry arguments one 
should start with the implicit constitutive law 






(e^^+j^iD-E + i^i{D-Ef + ..)E,, 

(' 



-il^B-H + a{B-Hf 



(31) 

(32) 



After this modification, all formulae of standard TO 
again can be used straightforwardly to arrive at the trans- 
formation medium 



D' = 



B' 






V7 



X^D-E- 



&BH 



^Z'iD-Ef 






renB-H) 



Ej. 

(33) 
H,. 

(34) 



These equations represent the simpHfied version of Eqs. 



(21) and (22 1 in the case where the nonhnear response 



can be expanded in terms of the energy densities rather 
than the electric and magnetic field monomials. 



It should be realized that the models (23l/(24) and 



(31)/(32) in general are two different limits of the the 
generic expansion (15)/(16). To illustrate this fact let us 
assume a specific situation where cubic and higher non- 
linearities can be neglected. Clearly, the model (23|/(24) 



corresponds to the situation, where the third and fourth 



order tensors in (15)/(16l may be written as 



ijk 

Xi 



il 



ijk 







ijkl 

Xi 



xi^ 



j-ijkl 
?2 



ij^,kl 



^2 7 



Conversely, (15)/(16l correspond to the limit 

ijk /-ijk n ijkl ^ij kl /-ijkl /^ij k 



(35) 



(36) 



Clearly the two descriptions are equivalent if and only 
if 7*-' ex £*^ ex /i*-', in particular in the case of the Kerr 
nonlinearity discussed below. 



B. The Kerr nonlinearity 

The theory in the previous subsection encompasses 
large classes of nonlinear media. In this subsection a fur- 
ther restriction to the special case of the Kerr nonlinear- 
ity is made. Kerr media are described by the constitutive 



D' = e 



XE') Y^E, , 



B' = ^ioY'H, 



(37) 



which, according to the discussion in the previous section, 
should be rewritten as 

D' = {e + xD- E) f^E, , B' = ^ial''H, , (38) 

with 

^ X ■ 

X = — h higher-order contributions, (39) 

and it is assumed that quartic and higher-order contribu- 
tions to the nonlinearity term are irrelevant. Taking the 
standard Kerr medium as a starting point, a coordinate 
transformation x ^ x yields a transformation medium 



D' = 



77 J V7 



^f'H, 



(40) 



Besides the transformation of the nonlinearity itself it 
is important to realize that the resulting medium has a 
non-trivial, though linear, magnetic response. 

Instead of the constitutive law (37 1, implications of 



the optical Kerr effect are often discussed in terms of an 
intensity dependent refractive index profile 



n-2 



no + ^2/ 
12^ 

2 a! 5 



no 



Stt 



(41) 
(42) 



where £ is the Fourier component of a monochro- 
matic wave E{t). In terms of this relation several self- 
interaction effects can be described, see e.g. Ref. [H]. As 
a specific example, let us briefly summarize the effect 
of self-focusing, whose behavior in TO is studied below. 
For powers of the incident beam that are greater than 
the critical power 



P., 



0.3727rAg 
8non2 



(43) 



with Aq being the vacuum wave length, self-focusing oc- 
curs with the self-focusing angle and self-focusing dis- 
tance (ignoring diffraction effects) 



2t12/ 

no 



Wq 



(44) 



Here, wq is the radius of the incident beam. If the beam 
power is equal to the critical power, the balance of self- 
focusing and diffraction effects results in self-trapping of 
the beam. It is important to note that self-focusing or 
self-trapping are specific consequences of the solution of 
the Maxwell equations in the Kerr medium (37). Since 



this solution maps algebraically onto the solution in the 
transformation medium ( 40 1 , the same effects will be 



seen in this new medium. In fact, by a suitable choice 
of transformations, self-interaction effects can be pre- 
cisely engineered using the TO approach. At the same 
time it should be realized that the description of self- 
interaction effects based on intensity dependent refrac- 
tive index profiles (such as Eqs. (43) and (44)) are not 



covariant and thus in general do not hold in the transfor- 
mation medium. Indeed, while in the Kerr medium only 
a scalar linear and nonlinear response are present, the 
more complicated transformation medium ( 40 ) in general 
is based on anisotropic material parameters. Quantities 
like a scalar refraction index thus are not defined in these 
media. 



EXAMPLES OF TRANSFORMED KERR 
MEDIA 



^ Evidently, in the special case of the Kerr medium the electric 
energy density is proportional to £^, which justifies the standard 
approach. However, for anisotropic media this ansatz remains 
problematic. 



Until now we have mainly developed the generic for- 
malism of nonlinear TO. In this section we discuss some 
explicit examples. As mentioned in the introduction al- 
ready, two different fields of applications of nonlinear TO 
exist: 



Nonlinear TO can be used in exactly the same way 
as linear TO with the sole difference that the start- 
ing point in this case is a nonlinear medium. An 
example of this type is discussed in Section |V A[ 

Nonlinear media exhibit interesting self-interaction 
effects, which can be engineered by means of non- 
linear TO. In contrast to the previous application, 
these devices rely on a productive combination of 
the self-interaction effects and the imposed coordi- 



With this result the relative permittivity and permeabil- 



nate transformation (see Section V B ) 



A. The nonlinear field concentrator 

As mentioned in the introduction, many different ap- 
plications of linear TO have been proposed. Most of 
these applications also work in nonlinear TO, with the 
important difference that in nonlinear TO the solutions 
in a Kerr medium are transformed as opposed to the 
free space solutions in standard linear TO. As a specific 
example we consider the field concentrator, whose trans- 
formations in cylindrical coordinates (r, </>, z) are [5] 



f r , < r < i?2 ; 




k^^ k^R^^ i?2<r<i?3; 


(45) 


r , r>R3; 




4> = 4> ; z ^ z ; 


(46) 



with < i?i < i?2 < -Rs- This transformation com- 
presses a cylinder with radius R2 in the original space 
into a cylinder with a smaller radius -Ri , conversely space 
is stretched in the region Ri < f < R^. In this way, 
the electromagnetic fields are concentrated in the cen- 
tral cylinder with radius i?i . Applying these transforma- 
tion rules to the inverse metric in cylindrical coordinates, 
[7'-'] = diag(l, l/r^, 1), one obtains 



ity follow from Eq. (11 1 as [5] 



[eH] = [f^] 



diag (1, 1, if j , < r < i?i ; 

g'^dmg(j„^,l) , R^<r<R,■, (50) 
^diag(l,;^,l) , r>R: 



3 • 



Here, we have used the abbreviation g = (i?3 — i?2)/(i?3 — 
i?l). From this result the nonlinear field concentrator 
follows immediately. From Eq. ( 38 ) the linear terms 
are given by the expression ( 50 1 - up to a rescaling of 



permittivity by the relative permittivity of the original 
medium in Eq. ( 38 ) . The transformation of the nonlinear 



term is slightly different, since this part is not rescaled 
by the determinants of the metrics. Thus, as the only 
missing piece 






diag(l,4,,l) , 



< r < i?i ; 
i?i < r < i?3 ; 
r >i?3 ; 

(51) 
is found as the relevant equation of the nonlinear field 
concentrator. 

The nonlinear field concentrator compresses the elec- 
tromagnetic field lines in a cylinder with radius R2 into 
a cylinder with smaller radius i?i, conversely fields are 
stretched around this smaller cylinder. This working 
principle is completely independent of the specific elec- 
tromagnetic fields in the material. 



B. Engineering self- focusing devices 



[f. 



where 



diag(^(^;^_^^)2,^(-)2, 
.diag(l,Jj,l), 



i-\ t-\ -R3 ^ -R2 . 
p(r) — r[r) ~ -^ ^^r 



0<f<Ri; 
i?i < f < i?3 ; 
f>i?3 ; 



i?2 — -Rl 



i?3 



R3 — Ri i?3 — i?i 

The determinants of the metrics are found as 

-ftr, 



Vl = r , 



R- 



-r 



^ ^ R3-R2 
R3-R1 
.f , 



< f < i?i ; 
p(f) , Ri <f < R3 

f>R3. 



(47) 
(48) 



(49) 



The above example does not make use of the special 
behavior of nonlinear media. Similarly to linear TO de- 
vices, its main feature (namely, concentrating fields in 
the cylinder with radius i?i) applies independently of 
the specific form of the electromagnetic fields which are 
transformed. However, an important aspect of nonlin- 
ear media is their interesting response to specific elec- 
tromagnetic field configurations. An interesting question 
thus arises, whether such effects also can be engineered 
using TO. In the second example we illustrate the power 
of TO in such an application as a self-focusing device. 
To keep the medium parameters at an acceptable level 
of complication, we require all medium tensors to be di- 
agonal matrices in Cartesian coordinates, but still allow 
them to be anisotropic. For a pulse propagating in the z- 
direction an interesting class of transformations obeying 



xiQ ^[m] 




E [V/m] 



xio ^[m] 



xio '[m] 




E [V/m] 



I 



I 



XIO ' [m] 



FIG. 1. Illustration of the laser beam in free space (A) and in a Kerr medium without transformation (B). Vertical lines indicate 
the region where the transformation will be applied in the second step, in this plot they are presented for illustrative purposes 
only. 



this constraint ifQ 



a 



-z , 



a > . 



(52) 



If the original coordinate system (x, y, z) is assumed to 
be Cartesian, then the above transformation results in 
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/7 = a , 



(53) 



and thus we find the transformed medium parameters as 
D' = (ae + xD-E\ f^Ej , B' = afiof^ Hj . (54) 



The parameter a defines how the typical length of self- 
interaction effects (e.g. the self-focusing length) is short- 
ened {a > 1) or streched (a < 1) with the transforma- 
tion. If we wish to express the nonlinearity purely in 
terms of Ei, we can make use of the relation 



oo 

D E = D'E, = aeY, T'^E^f^Ej] 

n=l 



E ^^'"-^ 



n=l 



(El + El 



'^ElT 



(55) 



Next, using a time-domain simulation in COMSOL 
Multiphysics [22], the possibility of engineering the self- 
focusing distance is demonstrated numerically. The re- 
sults are shown in Figs. [l}|3] Some explanations are in 



* More complicated manipulations are possible at the price of an 
inhomogeneous, nonlinear medium with the transformation 

^ = fiz), (f^) = diag(l,l,(/')2) , v^ = l//', 

where the function f(z) should be monotonically increasing. 



order: The left boundary of the figures is the excitation 
boundary, where a "laser beam" is modeled by the Gaus- 
sian beam 



y 



E{t) = £;oCxp ( -|^ ) sin(27ri/i) , 



(56) 



with the choice i?o = 1 V/m, ^ = 3 • 10^ ^'^ m and v = 
300 THz. The electric field points in the x-direction and 
propagates in the positive z-direction, i.e. from left to 
right in all figures. The side length of the plotted area is 
1.5 • 10"^ m in y and 2 • 10^^ m in the z-direction, with 
the snapshot taken at t = 10"^'' s; these numbers have 
been chosen for convenience and do not exhibit a deeper 
physical interpretation. The propagation of the beam in 
free space and in a Kerr medium 






(eo + XE') f'E, , 



X = 1.5 • 10" 



^10 



Cm 
V^ 



(57) 
(58) 



is illustrated in Fig. [T] This Kerr medium is used as a 
starting point of our manipulations; in the central re- 
gion of the simulation area, indicated by the two vertical 
lines, transformations with various values of a are cho- 
sen, while in the left and right areas all medium param- 
eters are kept fixed. Results for a — 0.25 and a — 2 are 
shown in Fig. [2] Since COMSOL Multiphysics does not 
allow for implicit definitions of the constitutive equation, 
the expansion ([55]) was used for the Kerr nonlinearity, 
whereby all the terms except for the leading one have 
been dropped. Despite this approximation, it is clearly 
seen that the self-focusing length is shortened for a > 1 
and elongated for a < 1. Furthermore, the interfaces be- 
tween the untransformed and the transformed Kerr me- 
dia are reflectionless; it can be shown on general grounds 
[2? that this must be true for any transformation of type 



(52 1 at interfaces defined by a constant coordinate z. If 
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FIG. 2. The effect of transformations with a = 2 (A) and a — 0.25 (B). The transformation is apphed to the central region 
indicated by the two vertical lines only, the left and right parts remain unchanged compared to the Kerr medium in Fig. IT] 



the self-focusing point lies in the transformed region (as 
is the case for a = 2), the self- focusing pattern is de- 
formed by the anisotropy of the medium parameters ( 54 ) 



(space is squeezed in the z-direction but not in the x and 
y-directions) . 

Completely new effects can be obtained by choos- 
ing extreme medium parameters, for example a very 
small value of a can yield a freezing of the focal points. 
To demonstrate this we change the untransformed Kerr 
medium in such a way that the focal point is located ex- 
actly at the interface of the left and center areas (see Fig. 
^. In the second step a transformation is again applied 
to the central region, in this example with a very small 
value of a = 0.01. In this way the width of the cen- 
tral region becomes very small in the transformed space, 
the self-interaction effects can hardly evolve and conse- 
quently the self-focusing is frozen. As a result we obtain 
an almost ideal flux tube over the whole area of the trans- 
formation medium. 

This example shows that genuine characteristics of 
nonlinear media can be engineered by means of nonlinear 
TO. Of course, such devices will only work for a class of 
electromagnetic fields that exhibit a certain characteris- 
tic, in contrast to the examples of the previous subsec- 
tion, whose workability is completely general. Still, this 
example shows a completely new field of applications for 
nonlinear TO, which does not exist in conventional linear 
TO. 



VI. CONCLUSIONS 

In this paper we have developed the formalism of non- 
linear transformation optics. Standard transformation 
media, which can be interpreted as mimicking empty 
space, are inherently linear. Therefore, in nonlinear TO 
the original space, on which the transformation is ap- 
plied, has to be nonlinear. In a generic expansion of 



the constitutive relations in terms of E and H the for- 
malism extends straightforwardly from the linear case, 
at least for purely spatial transformations. Less imme- 
diate are situations such as the Kerr nonlinearity, where 
the nonlinear contributions are expanded in terms of E'^ 
and H^, since these quantities implicitly depend on the 
metric: E^ = Eij^^Ej. We have shown that this prob- 
lem can be circumvented by defining implicit constitutive 
relations in terms the electric and magnetic energy den- 
sities, D ■ E and B ■ H. In this formulation nonlinear 
TO retains all the intuitive design properties of standard 
linear TO, such that the (nonlinear) field configurations 
simply can be considered as sticked to the coordinate 
lines during the transformation. 

Most applications of linear TO are insensitive to the 
exact form of the electromagnetic fields, in other words, 
their functional principle works - at least theoretically - 
for any field configuration. Such devices extend straight- 
forwardly to nonlinear TO, which was demonstrated in 
an example of an electromagnetic concentrator in a Kerr 
medium. Besides such applications, nonlinear TO allows 
for a completely new type of devices: Since nonlinear 
media exhibit interesting self-interaction effects, which 
depend on the exact field configuration, TO can be used 
to engineer such effects. As a specific example we have 
shown that the self-focusing effect of a laser beam in a 
Kerr medium can be engineered by means of nonlinear 
TO. Although the working principle of such devices is 
restricted to exactly those field configurations that ex- 
hibit the desired self-interaction effect in the original Kerr 
medium, the design principles of TO are retained: By 
a suitable compression or stretching of space the self- 
focusing length can be altered; if space is stretched by 
a large factor exactly at the focusing point, the focusing 
area is "frozen" , which results in a high intensity flux 
tube. 

The phenomenology of nonlinear media is much richer 
than we could treat in this article. It would be interesting 
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FIG. 3. The untransformed Kerr medium with x = 1-9 ' 10 ^ Vfi^ (A)- On the right hand side (B), a — 0.01 is chosen in the 
central region indicated by the two vertical lines, such that the focal point freezes and builds a concentrated flux tube. 



to study the implications of nonlinear TO in more detail, 
including such aspects as the global solutions (solitons), 
second harmonic generation or different self-interaction 
effects. On a different route of research it should be in- 
vestigated how the proposed nonlinear materials could 
be realized as artificial media (metamaterials) . 
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